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ABSTRACT 

We present HST/STIS time-series spectroscopy of the central star of the 'Cat's Eye' planetary nebula 
NGC 6543. Intensive monitoring of the UV lines over a 5.8 hour period reveals well defined details of large- 
scale structure in the fast wind, which are exploited to provide new constraints on the rotation rate of the central 
star. We derive characteristics of the line profile variability that support a physical origin due to co-rotating in- 
teraction regions (CIRs) that are rooted at the stellar surface. The recurrence time of the observed spectral 
signatures of the CIRs is used to estimate the rotation period of the central star and, adopting a radius between 
0.3 and 0.6 Rq constrains the rotational velocity to the range 54 < Vrot < 108 kms~^ The implications of these 
results for single star evolution are discussed based on models calculated here for low-mass stars. Our models 
predict a sub-surface convective layer in NGC 6543 which we argue to be causally connected to the occurrence 
of structure in the fast wind. 

Subject headings: stars: winds, outflows — stars: evolution — ultraviolet: stars 



1. INTRODUCTION 

Planetary nebulae (PNe) display remarkably diverse mor- 
phologies, the large majority of which are not spherical. The 
formation and shaping of bipolar PNe, which can include axi- 
and point- symmetries, is currently a moot issue in stellar as- 
trophysics, with the likelihood that different mechanisms may 
be more dominant in different nebulae. An understanding 
of this PNe diversity encompasses the potential roles of sev- 
eral interesting physical processes, including magnetic fields, 
common envelope evolution in binary systems, stellar radia- 
tion fields and outflows. 

One of contending scenarios for sculpting the nebulae is 
based on stellar rotation in the Generalized Interacting Stellar 
Wind model (GISW e.g., Frank 1999). In this case, asymme- 
try results from an equatorially enhanced outflow from a post- 
asymptotic giant branch (AGB). The higher equatorial density 
constrains the nebular expansion in that region, leading to an 
elliptical or bipolar PN. The origin of the requisite equatorial 
expansion remains uncertain, but rotationally induced aspher- 
ical mass-loss from the central star is one of the favoure d sce- 
narios (e.g., Dwarkadas & Owocki 2002; Garcia-Segura et al] 
[2005). In a radiatively driven wind, rotation may result in lat- 
itudinal variations in the wind, with high velocities and larger 
wind-momentum gas at the poles compared to the equator. 
Rotation can also significantly affect the magnetic channelling 
of a hot star wind (e.g., Ud-Doula et al. 2008). 

A pivotal issue in the role of stellar rotation in PNe - 
and the subject of this Letter - is whether the precursor 
(AGB) and remnant (CSPN) rotation is fast enough to pro- 
duce the contrasting equatorial and polar outflows. Observa- 
tional constraints on AGB rotation rates suggest upper limits 
of a few km s~^ (e.g., de Medeiros & Ma vor 1999), though 
llgnace et al.l fl996.) argue that very modest rotation rates in 
AGB stars are not likely to result in aspherical geometries. 
In their study of the ev olution of r otating stars from main- 
sequence to white dwarf. ISuijs et al.l ([200 8) predict that single 
star progenitors of CSPNe have very little or no rotation at all 



in this phase. 

In this Letter we provide new perspectives and constraints 
on the rotation of CSPN by studying the temporal evolu- 
tion of large-scale (coherent) wind structure, which we ar- 
gue is rooted at the stellar surface. Over the past two 
decades extensive UV spectral time-series observations of 
OB and WR stars, obtained primarily with the Interna- 
tional Ultraviolet Explorer (lUE) satellite, have revealed 
systematic behaviour of wind variability. The most thor- 
oughly studied objects exhibit cyclical variation in the stel- 
lar winds that are causally connected to the photosphere, 
with the time- scales of the variability r elated to the rota- 
tion periods of the mas sive star s (e.g., 'Massa et al.l 1 1995 : 
Fullerton et al. 1996; Ka oeretaLl 1996: deJongetal. 200TI: 
.Prinja et al.. .2002: .Chene & St-Lo uis 2011). The physical 
interpretation of this rotational modulation is the action of 
spiral- shaped structure in the wind flow caused by perturba- 
tions at the base of the wind and carried by stellar rotation. 
Using Far- Ultraviolet Spectroscopic Explorer (FUSE) data, 
similar wind structure ha s been demonstra t ed to exist in the 
fast winds of CSPN fe.g. JPrinja e'tani2007L[20T2l) . 

In this Letter we report on high-resolution, high-quality 
Hubble Space Telescope (HST) time-series spectroscopy of 
the central star of the hydrogen-rich (Cat's Eye) nebula 
NGC 6543. The time dependence of the fast wind of 
NGC 6543 was previously shown to exhibit discrete absorp- 
tion components, which traverse the absorption trough s of 
unsaturated wind lines on time-scales < 1 hour (Prinja e t al.l 
,2007). We present here a multi-line study of the variability in 
the HST spectra. The clarity of the structure seen in the res- 
onance, excited state and iron UV lines provides key details 
on the nature of the wind structure, and thus its origins and 
connection to the rotation of the central star. 

2. //OTSTIS OBSERVATIONS 

The central star (HD 164963) of the well studied PN 
NGC 6543 was observed with HST for Program #12489 (D. 
Massa PI) on 2012 May 16. Repeated Space Telescope Imag- 
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TABLE 1 

Adopted parameters for the central star of NGC 6543 



Parameter 


Value 


Spectral type 


Of(H-rich) 


Luminosity 


1585L0 


^eff 


67000 K 


Radius 


O.3R0 


Radial velocity 


-66 km s~^ 


Terminal velocity 


1340kms-i 


Velocity law 


2 



Note. 
2008 



— Parameters are from the non-LTE analysis of Georgiev et al. 
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Fig. 1. — Dynamic spectra revealing evidence for systematic variability 
in almost every detectable absorption line in the HST/STIS spectrum of the 
central star of NGC 6543. Recurrent large-scale wind structure is diagnosed 
by the Siiv resonance doublet and Oiv and Ov excited state lines. Sympa- 
thetic variations are also evident in numerous Fevi (AA1260 to 1310) and 
Fev (AA1390 to 1460) lines. 

ing Spectrograph (STIS) E140M grating (spectral resolution 
10 kms~^) observations were secured for 5 contiguous HST 
orbits in ACCUM mode and with individual exposure times of 
180-s. A total time-series of 95 spectra were secured spanning 
~ 5.8 hours, each with a wavelength range ~AA1140tol710 
A. The signal-to-noise resolution per element of an individual 
spectrum varies from ~ 5 at the extreme wavelength ends to 
-30at 1300A. 

The rich UV spectrum of HD 164963 contains several key 
diagnostic lines including the resonance line doublets of Nv 
AA1239, 1243, Siiv AA1394, 1403, and Civ AA1549,1551; 
the excited state wind lines of Oiv AA1339, 1344 and Ov 
A1371; Hell A1640, as well as 'forests' of Fev and Fevi ab- 
sorption lines. Table 1 lists the fundamental parameters of 
HD 164963 we adopt for this paper. 

3. EVIDENCE FOR EVOLVING LARGE-SCALE WIND STRUCTURE 

Previous studies established that the fast wind of NGC 6543 
is variable on short (~ hourly) time- scales, with diagnos- 
tic ch anges evident in optical spec troscopy (IMendez et al.l 
[T99Qh . Dh otometrv feell et al.| [T994h and FUV to UV spec - 
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troscopy (iPrinja et al.l 
The widespread and rapid line profile variation present in our 
HST/STIS time- series is demonstrated in Fig. 1 which shows 



ratios of the individual STIS spectra divided by the mean for 
the time- series as a dynamic spectrum (i.e. a two-dimensional 
image of the time ordered ratios). 

The temporal behaviour seen in Fig. 1 is remarkable. In ad- 
dition to the 'classic' blueward migrating and recurring dis- 
crete absorption components (DACs) present in the unsatu- 
rated resonance line doublet of Siiv AA1394, 1403 and in 
the subordinate excited state lines of Oiv AA1339, 1344 and 
OvA1371, the numerous Fevi absorption features between ~ 
AA1260 to 1310A and Fev lines between AA1390 to 1460A 
also vary in concert with the wind lines. Three DAC episodes 
are evident over 5.8 hr (see also Fig. 2): There is some em- 
pirical resemblance in optical depth and velocity dispersion 
between the first episode (T ~ 0.3 to 2 hr) and the third fea- 
ture (T - 3.7 to 5.5 hr). The second DAC episode (T - 2.3 
to 3.3 hr) is weaker though it has a similar blueward velocity 
progression. We note that the individual spectra were normal- 
ized to suppress an overall continuum variability of ~ ±4%, 
which is similar in amplitude and time scale to that seen in the 
optical by Bell et al. (1994). Because the time dependence is 
not obviously in phase with the DACs, its exact origin is not 
clear. However, the wavelength independence of the contin- 
uum variability suggests that it results from some sort of geo- 
metric obscuration instead of temperature irregularities. 

The reaction of the deep-seated Fe lines indicates that the 
wind structures originate very close to or at the stellar surface. 
The variations in excited lines of Oiv and Ov provide further 
strong support of this deep-seated origin. In general, a stel- 
lar flow need not be a monotonically increasing function of 
distance from the star and an absorption feature seen at low 
velocity i n a resonance line need not originate close to the 
star (e.g. JOwocki et al.lll988l) . In contrast, in dense expand- 
ing winds the populations of the excited levels require a strong 
EUV radiation field, thus these transitions have a strong de- 
pendence on the radial distance from the star. This means that 
large populations of the lower levels of Oiv and Ov ions can 
only occur close to the star. We conclude therefore that the 
evolving features seen in Fig. 1 in the fast wind NGC 6543 
originate at the stellar surface. 

The short time-scale fluctuations in the numerous Fe 'for- 
est' lines also imply that the surface conditions on the CSPN 
are non-uniform, and likely to affect the results of 1-D model 
atmosp here analysis . This may account for the inconsisten- 
cies Georgi ev et al.l (12008) report in their non-LTE analysis 
of HD 164963, including an overabundance of Si, observed 
spectral features not present in the model, differing wind ve- 
locities for different resonance lines, and the need for addi- 
tional ionization sources in the external part of the fast wind. 

4. CO-ROTATING INTERACTION REGIONS 

The variability patterns observed in OlV, Ov and 
Hell A 1641 are shown in closer (velocity) detail in Fig. 2, 
where some further characteristics are revealed. Firstly, the 
line profile modulations tend to be bow- shaped in the di- 
rection of increasing time, i.e., an absorption enhancement 
occurs first at an intermediate velocity and then spreads to 
higher (bluer) and lower (redward) velocities simultaneously 
at a later time. This 'phase bowing' is particularly evident in 
the dynamic spectra of OlV and Hell. Similar forms of vari- 
abilit y have been identified in the stellar winds of OB stars 
(e.g. iFullerton et al.l 119971: lie Jons et al.' '20011 ; iPrinia et al.l 
120021) . H ydrodynamic mode ls by Cranmer & O wockil ([l996l) 
(see also 'Lobe l & Blommd 12008^) showed that large scale, 
pro-azimuthal, wind structure in the form of co-rotating in- 
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Fig. 2. — Dynamic spectra of Oiv, Ov and Heii revealing bow-shaped 
modulation that extends blueward and redward simultaneously in time. Note 
also the ionization shift between Ov and the lower ions Oiv and Heii. 

teraction regions (CIRs) can explain the characteristics of the 
bowed discrete absorption components. The CIR structures 
are akin to those seen in the solar wind (e.g., iHundhausenI 
0^972; Mullan 1984), where different sectors of the stellar sur- 
face give rise to winds which accelerate differently. When, 
due to rotation, these different sectorial flows interact, an in- 
terface with a spiral pattern results and a velocity plateau de- 
velops along the interface. As the CIR pattern rotates across 
our line of sight to the stellar disk, it creates an absorption 
in the dynamic spectrum which evolves toward high veloc- 
ity. Intensity irregularities at the stellar surface, such as spots, 
non-radial pulsations, or granulation can potentially alter the 
radiative wind acceleration locally, thus giving rise to streams 
of faster and slower material, which interact to form the CIRs. 

We also note in Fig. 2 evidence for an ionization shift, 
where the absorption enhancements with respect to the mean 
profile (i.e. darker patches) in Ov lag those in the OlV and 
Hell by about 1 hr. (See, e.g., the episode between AT ~ 4 
and 5-hr in Fig. 2). This may be evidence for an ionization 
gradient across the CIR, with the lower ions on the leading 
edge of the CIR and higher ions on the trailing edge. 

5. CONSTRAINING THE STELLAR ROTATION RATE 

We have demonstrated that the UV wind lines of NGC 6543 
are extensively variable, and that the variability results from 
large-scale, coherent wind structure. The characteristics of 
the variability, including bow- shaped modulations, are consis- 
tent with co-rotating interaction regions (CIRs) that are rooted 
at the stellar surface. In this section we use the temporal evo- 
lution of the DACs to constrain the rotation rate of the central 
star in NGC 6543. 

By their very nature, CIRs infer rotation. They have been 
detec ted in t he winds o f O ( Kaper et al. 1997; de Jong et al] 
1200 Ih B (Fullerton et al. 1997rPrinja et aL| |2002|) and WR 
stars (e.g., Chene & St-Louis 2011). Typically, the signatures 
of 2 strong CIRs are seen, but these are sometimes accom- 
panied by two weaker arms (Fullerton et al. 1997, de Jong 
et al. 2001). It seems reasonable, therefore, to associate the 
two strong features which occur at c^03 and 3.7 hours 
in Figures 1 & 2 as strong CIRs and the one near 2.3 hours 
as a weaker one. Thus, we estimate the rotation period of 
HD 164963 is 2 x (3.7-0.3) - 6.8 hours. In this case, its 
rotation rate is ~ 1797? km s~\ where R is the stellar radius 
in solar units. The most recent line-blanketed non-LTE line 
synthesis analysis of the central star yields estimates of R for 



HD 164963 that vary between 0.3 (Table 1) and 0.6 (Herald 
and Bianchi 2011). We therefore estimate its rotation rate to 
be in the range 54 < Vrot < 108 kms~^ 

6. DISCUSSION 

State-of-the-art stellar evolution calculations of rotating, 
single low-mass stars predict a substantial spin-down dur- 
ing their late phases of evolution. This is because, despite 
a possibly high initial rotation rate, the dramatic expansion 
and increased mass-loss during the RGB and AGB phases 
result in huge angular momentum loss. Coupling between 
the core and the envelope is provided by angular momentum 
transport processes, arising from rotational instabilities and 
magnetic torques. This coupling is required to explain the 
slow rotation rate of wh ite dwarfs and, for more massive pro - 
genitors, neutron stars (iHeger et al.l 120051: ISuijs et al.ll2008h . 
Moreover an angular momentum transport in radiative zones 
like the one provided by magnetic torques is necessary to re- 
produce the near uniform ity of the solar rotation profile (e.g., 
lEggenberger et al.l 120051) and possibly the recent asteroseis- 
mic observations of the red g iant KIC 8366239 (iBeck et al.l 
1201 2l: lEggenber ger et al.|[2012 ). However we want to stress 
that the angular momentum transport produced by magnetic 
torques relies on the existence of a dynamo in the radiative 
layers of stars (Tayler- Spruit dynamo. S pruit 2002), which is 
currently debated (e.g. JZahn et al.|[2007l) . Other angular mo- 
mentum transport processes i n radiative zones have also been 
proposed (e.g., gravity waves Pl^on & Charbonnelll2005|) . 

We have used the Modules for Stellar Experiments in As- 
trophysics (Mesa version 4298 Paxton et al. 201 1) code to cal- 
culate the evolution of a grid of low-mass stars (1.5M0, 2Mq 
and 3Mq) with metallicity Z = 0.02 and initial equatorial rota- 
tional velocity of 250kms~^ . This code accounts for the trans- 
port of angular momentum by rotational instabilities and mag- 
n etic torques (see e.g., Hegeretal. 2000, 2005; Suiis etaD 
|2008|) . We adopted the mass-loss rate of Reimers (1975) dur- 
ing the RGB evolution and the one of Bloecker (1995) during 
the AGB and post- AGB phase. Our 1.5M0 model provides 
the best fit to the values of luminosity ( 1 5 8 5L(7^ ) and effectiv e 
temperature (67000K) determined bv lGeorgiev et al.1 (l2008h . 
however the major results discussed below do not change for 
the other calculations. We find that when the star reaches the 
CSPN phase, it has lost most of its angular momentum and 
its surface rotation velocity is much lower than lkms~^ (see 
Fig. 3). Therefore calculations that adopt the physical in- 
gredients required to explain both progenitor and descendant 
stars of CSPN, are at odds with the limits on the rotation rate 
of NGC 6543 derived in Sect. 5. Note tha t DACs seem to 
be frequent among CSPN (iPrinja et al .1120 12 ). Since to be ob- 
servable DACs usually need to recur on a short period, this 
suggests that rapid rotation could be widespread among CSPN 
(where here rapid means rotating faster than the predictions of 
single stellar evolution calculations). 

An interesting possibility is that NGC 6543 did not evolve 
as a single star. In a binary system spin-up of one or both 
components could result from mass accretion, tidal interac- 
tion or a merger. This is particularly interesting in the case 
of CSPN, as binarity has been often discussed as a possible 



explanation for the asymmetric shapes of PNe (see e.g., IbenI 
1995; Balic k & Franl200l: lluarte-Espinosa et al.ll2012D . 



CSPN are hot and luminous, and due to an iron opacity peak 
at about 150000 K, they can develop a subsurface convection 
region. This occupies a negligible fraction of the stellar mass, 
but at the same time occurs very close to the stellar surface. 
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Fig. 3. — Evolution of the equatorial rotational velocity for our 1.5 Mq 
model with metallicity Z = 0.02 and initial equatorial rotational velocity of 
250kms~^. The horizontal log scale shows the time to the end of the calcu- 
lation, when the luminosity and surface temperature roughly correspond to 
the estimated values for the central star of NGC 6543 (See Table 1). The 
model includes angular momentum transport induced by rotational instabil- 
ities and magnetic torques. The evolution of stellar luminosity (dotted line) 
and surface temperature (dashed line) are also shown for reference. 

Our model of the central star of NGC 6543 reveals the pres- 
ence of such convective layer, with convective velocities of 
about lOkms"^ (~15% of the local sound speed; CSPN mod- 
els with higher initial mass result in higher velocities). Note 
that the presence of this convective layer is a robust result. 
It is a structural property depending only on the opacity and 
luminosity, and therefore we do not expect major changes if 
the star went through binary interactions. Like the case of 



massive OB stars, it is possible that such convection zones 
can produce density and velocity fluctuations at the surface, 
stochastically excite oscillations and pr oduce magnetic fields 
(ICantiello et al.ll2009l: ICantiello & Bra ithwaite 2011). Inter- 
estingly the recent observation of stochastic brightness vari- 
ations in the CSPN of NGC 6826 ( Jevtic et al. 2012) could 
be related to the presence of such subsurface convective lay- 
ers. The non-uniform surface conditions revealed by the FeV 
and Fe VI line profile variabilityin NGC 6543 seem to further 
support this notion. If a dynamq^is at work in the iron convec- 
tion zone of NGC 6543, adopting the densities and MLT con- 
vection velocities from our models one finds an equipartition 
field B~2kG (inside the convective region). Similar to the 
case o f massive OB stars discussed in lCantiello & Braithwaitd 
(l2QTTh . these fields can in principle reach the stellar surface 
through magnetic buoyancy and affect the stellar surface (e.g. 
creating bright spots). Bright spo ts have been discussed as 
a way to produce CIRs (C ranmer & Owockil [1996). which 
might provide a n explanation for the co mmon emergence of 
DACs in CSPN dPrinja et al.ll2"QQ7[l2QT2h . 

Our work shows that observations of DACs in CSPN can 
provide important constraints on the surface rotational veloc- 
ity of these intriguing stars. In turn this information can be 
used to constrain their evolutionary history, potentially shed- 
ding light on the puzzling origin and shaping of planetary neb- 
ulae. 

MC would like to thank Bill Paxton for useful discussions. 
Support for this work was provided by NASA through grant 
number HST-GOO 12489 from the Space Telescope Science 
Institute, which is operated by AURA, Inc., under NASA con- 
tract NAS5-26555. This research was supported in part by the 
National Science Foundation under Grant No. NSF PHY05- 
51164. 



REFERENCES 



Balick, B., & Frank, A. 2002, ARA&A, 40, 439 

Beck, P. G., Montalban, J., Kallinger, T., et al. 2012, Nature, 481, 55 

Bell, S. A., Pollacco, D. L., & Hilditch, R. W. 1994, MNRAS, 267, 1053 

Bloecker, T. 1995, A&A, 297, 727 

Cantiello, M., & Braithwaite, J. 2011, A&A, 534, A140 

Cantiello, M., Langer, N., Brott, I., et al. 2009, A&A, 499, 279 

Chene, A.-N., & St-Louis, N. 2011, ApJ, 736, 140 

Cranmer, S. R., & Owocki, S. R 1996, ApJ, 462, 469 

de Jong, J. A., Henrichs, H. F, Kaper, L., et al. 2001, A&A, 368, 601 

de Medeiros, J. R., & Mayor, M. 1999, A&AS, 139, 433 

Dwarkadas, V. V., & Owocki, S. R 2002, ApJ, 581, 1337 

Eggenberger, R, Maeder, A., & Meynet, G. 2005, A&A, 440, L9 

Eggenberger, P., Montalban, J., & Miglio, A. 2012, A&A, 544, L4 

Frank, A. 1999, New Astronomy Review, 43, 31 

FuUerton, A. W., Gies, D. R., & Bolton, C. T. 1996, ApJS, 103, 475 

FuUerton, A. W., Massa, D. L., Prinja, R. K., Owocki, S. P., & Cranmer, 

S. R. 1997, A&A, 327, 699 
Garcia-Segura, G., Lopez, J. A., & Franco, J. 2005, ApJ, 618, 919 
Georgiev, L. N., Peimbert, M., Hillier, D. J., et al. 2008, ApJ, 681, 333 
Heger, A., Langer, N., & Woosley, S. E. 2000, ApJ, 528, 368 
Heger, A., Woosley, S. E., & Spruit, H. C. 2005, ApJ, 626, 350 
Huarte-Espinosa, M., Frank, A., Balick, B., et al. 2012, MNRAS, 424, 2055 
Hundhausen, A. J. 1972, Coronal Expansion and Solar Wind 
Iben, Jr., I. 1995, Phys. Rep., 250, 2 



Ignace, R., Cassinelli, J. R, & Bjorkman, J. E. 1996, ApJ, 459, 671 

Jevtic, N., Stine, R, Nilsen, W., et al. 2012, ApJ, 756, 9 

Kaper, L., Henrichs, H. R, Nichols, J. S., et al. 1996, A&AS, 116, 257 

Kaper, L., Henrichs, H. R, FuUerton, A. W., et al. 1997, A&A, 327, 281 

Lobel, A., & Blomme, R. 2008, ApJ, 678, 408 

Massa, D., FuUerton, A. W., Nichols, J. S., et al. 1995, ApJ, 452, L53 

Mendez, R. H., Herrero, A., & Manchado, A. 1990, A&A, 229, 152 

MuUan, D. J. 1984, ApJ, 283, 303 

Owocki, S. R, Castor, J. I., & Rybicki, G. B. 1988, ApJ, 335, 914 
Patriarchi, R, & Perinotto, M. 1997, A&AS, 126, 385 
Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3 
Prinja, R. K., Hodges, S. E., Massa, D. L., FuUerton, A. W., & Burnley, 

A. W. 2007, MNRAS, 382, 299 
Prinja, R. K., Massa, D., & FuUerton, A. W. 2002, A&A, 388, 587 
Prinja, R. K., Massa, D. L., Urbaneja, M. A., & Kudritzki, R.-P. 2012, 

MNRAS, 422, 3142 
Reimers, D. 1975, Memoires of the Societe Royale des Sciences de Liege, 8, 

369 

Spruit, H. C. 2002, A&A, 381, 923 

Suijs, M. R L., Langer, N., Poelarends, A.-J., et al. 2008, A&A, 481, L87 
Talon, S., & Charbonnel, C. 2005, A&A, 440, 981 

Ud-Doula, A., Owocki, S. R, & Townsend, R. H. D. 2008, MNRAS, 385, 97 
Zahn, J.-R, Brun, A. S., & Mathis, S. 2007, A&A, 474, 145 



1 



The convective turnover timescale is of order hours, comparable to the • ^^j^^ .Unt nr. ^, , ^,,r,nrv.^ rv.4„v,t nt 
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mf erred rotational period of NGC 6543. This means that the Rossby number 



